function Uout = multi2Q(Uin,np,emax,fieldeps,E0,muop)

%        multiplication step in the Chebychev reccursion
        
  
        Uout = hamil2Q(Uin,np,fieldeps,E0,muop);

%        normalization of the Hamiltonian operator range
%        from -1 to 1 :
%        H' =  ( 2 H - DE ) /DE 

        Uout = 2*Uout/emax - Uin ;
       end